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SUMMARY 


The present treatment investigates the transient behavior of an 
airfoil in supersonic’ flow due to pitching, flapping, and vertical gusts 
by means of the Fourier integral method. It is an extension of the 
earlier work "The Transient Reaction of an Airfoil Due to Change in 
Angle of Attack at Supersonic Speed" by Chieh-Chien Chang, Journal of 
the Aeronautical Sciences, volume 15, number 11, November 19*t8, 
pages 635-655* If a control surface is attached to the airfoil and 
deflected according to a prescribed time schedule, the effect can be 
evaluated from the result found for the main wing. An example is shown. 

By using a convolution integral as shown in the above-mentioned 
article, the aerodynamic effect of any arbitrary motion of the airfoil 
can be determined, if the aerodynamic response of the airfoil to a step- 
function disturbance is known. As useful examples, the harmonically 
oscillating airfoil with different modes or degrees of freedom is 
analyzed for the complete time history, for the case in which the motion 
starts abruptly from rest. The results check exactly with those of 
Garrick and Rubinow in NACA TN 1158, after the transient effect of the 

abrupt beginning has died out, that is, for times later than 

from the beginning of harmonic oscillation where M is the Mach number, 

U is the free-stream velocity, and C is the airfoil chord in con- 
sistent units. 

In carrying out the investigation, a new integral or function 
C(3,M) has been introduced. This new function is sufficiently funda- 
mental to the supersonic flutter problem that the behavior of all the 
different degrees of freedom can be expressed in terms of C(p,M) and 
the usual transcendental function (here (3 is related to the frequency 
of oscillation). It is believed that this new function plays an equiva- 
lent part in supersonic flutter to the Theodorsen function in the flutter 
problems of incompressible flow. Preliminary investigation of both the 
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complete C-function C(|3,M) and the incomplete 0 C -function is made. 

The incomplete C-function together with the incomplete Bessel function, 
so designated by Garrick and Rubinow, are important to the transient 

region for times between (—7— )£ and A n extensive calcu- 

\M + 1/U , \M - 1/U 

lation of this new function will be very useful in supersonic flutter 
analysis. 


INTRODUCTION 


Recently, the transient problem of the nonstationary motion of an 
airfoil in supersonic flow has been treated by a number of investigators. 
With the Fourier integral method the present author first presented a 
paper on the transient reaction of an airfoil due to change in angle 
of attack at supersonic speed (reference 1) in January 19^. As a 
byproduct of the investigation, the transient effect due to a vertically 
descending airfoil is also shown in reference 1 . Concurrently, Biot 
(reference 2) investigated Schwarz's problem (reference 3 ) of vertical 
gust with a different approach. Later Heaslet and Lomax (reference L) 
considered the cases of the vertically descending airfoil and an air- 
foil encountering a vertical gust with the frame of reference at rest 
with the undisturbed fluid. Miles treated the above gust and flapping 
case with Green's function (reference 5). 

The present paper is an extension of the author' s earlier work 
(reference 1 ) to all possible modes of nonstationary motion, and it 
shows that the results can be applied to the transient effect of con- 
trol surfaces. With the concept of the convolution integrals as 
developed in reference 1 , the transient results can be applied to many 
types of nonstationary motion, particularly the harmonic oscillation. 

There are four types of problems to be treated in this paper. 

(a) Investigation of the case of a suddenly descending airfoil 
with a constant rate - so-called flapping degree of freedom of the wing 

(b) Investigation of the case of the airfoil encountering a vertical 
gust of constant strength 

(c) Application of the convolution integral to evaluate the wing 
loading due to different modes of nonstationary motion, such as harmonic 
oscillations, and so forth 

(d) Extension of the present results to other degrees of freedom, 
such as aileron deflection, and so forth 
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The investigation of drag-characteristic behavior of the wing in 
steady flow is omitted entirely from this paper, because it is not a 
difficult matter for a designer to calculate the drag characteristics 
with the information contained in this paper and available solutions 
of the steady case. 

The present investigation was conducted at the Johns Hopkins 
University under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. The author would like 
to express his appreciation to Misses Vivian O'Brien and Patricia 
darken for their assistance in carrying out the project. 


' TRANSIENT REACTION OF AN AIRFOIL DUE TO UNIT-STEP CHANGE 
OF FLAPPING AND VERTICAL GUST 


•As shown in references 1 and 6, the linearized partial differential 
equation of the irrotational flow of a compressible nonviscous fluid is 



where 0 is the disturbance potential of the flow, the x-axis is 
positive in the opposite direction of the flight, and the y-axis is 
positive in the upward direction. The origin is attached to the air- 
foil leading edge. The details of the notations are given in 
appendix A. In the supersonic case, the disturbance potential of such 
an airfoil is equivalent to that due to a source sheet at the y = 0 
plane with time -dependent strength and is 


0U,y,t) 


1 1 f 11 | 2 v(l,+0,t - t) 

* (m 2 - l) 1 / 2 Jo Jt! ^|(t - r x ) (t 2 - t) 


( 2 ) 


where v(£,+0,t - t) is the vertical velocity component on the wing 
surface- at the point (£,+0) at an earlier time t - r. The term 

v(i,+0,t - t) . ■ ■ 

— is just the source-sheet strength per unit area. In 

( m 2 - i ) l/2 

addition, 
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If the point of interest is on the lower side of the wing surface 
(x,-0,t) .equation (2) reduces to 


0(x,-O,t) 



(5) 


With the above equation the transient case due to. sudden change of 
angle of attack of the airfoil about any arbitrary point has been treated 
in reference 1. The important results of that paper are given in table 1 
for use in the present development. ,Two additional interesting cases 
will be treated as follows. 


Vertical Flapping 

/ 

At time tQ, the airfoil suddenly descends with constant velocity fi. 
The unit-step downwash at any point on the wing can be expressed in terms 
of Dirichlet's integral as 

v(x,-0,t) = fi(t) = fi - f Sin 60 — — + 

_ n Jo ' 01 
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where o> is the angular velocity per second or a frequency parameter. 
Figure 1(a) shows ^ - = a^(t). Substituting equation (6) into equa- 
tion (5), there results, . after manipulation, 


IzfhU,- 0 ^) = 


* 2 (m 2 - l) 1 / 2 


PTo P 00 / \ 

. ■ sin 0) ft - r - to) 

d| dT — V ... . — — d£o + 

J T 1 Jo " T l )( T 2 ~ t ) 


Ox 


21 

2 


PTr 


Jo U 


dT 


T 1 ^( T - T l)( T 2 - T ) 


lx j'5 

h I I da> 


"Jo - sin “ (t ' t0 ’ " W '‘ ) 


hx 


2(m 2 - l) 1 / 2 


(7) 


or after integration the potential is expressed in three zones according 


to the value of. 


U(t - t 0 ) 


= v as follows: 


0 h (x,-O,t) = 


fix 

5T ' 


■ — - — - j- sin -1 f- - m) - 

(m 2 - 1)1/2 U ) 

■Lfsin- 1 M(1 - v) - *1 S V < 

Mv L 2j V M - - M / 


-it 


( m 2 - l ) 1 / 2 


( 


0 £ v ^ 


M - 
M 




TC 

Mv 


/ M + 1 < \ 
[~ S V 


( 8 ) 
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where 


v = 


U 


A 

fa " t o) (y - t o') 


* A 

X 1 = — 


m = 


M 2 - 1 


M 


f -S2 


X = 


i(m 2 - l) 


V = 


t 0 u 


and J 0 (cnX.) is the Bessel function of the first kind of the zero order 

K } it ■ 

“ a h 


Table 2 gives the detailed characteristics of — , . -E* 1 = f, ^h _ t? 

*r. ’ — ri’ 


Cm^ 

and ~ — = Gjj as functions of time. Figures 2 to 4 show these charac- 
a h 

teristics. The detailed methods of calculation are quite similar to 
those given in reference 1, and are omitted in this paper. 


In figures 2 to 4, a few important features may be explained. Within 
the zone III 

to Mach number and independent of time and the location of x' . But in 
time zone I 


•es 2 to 4, a few important features may be explained. Withii 

Ir, < t> " V < M \ „ C ph 

m — -s ——I, f h = - — is inversely proportional 

\ x M + 1/ a h 

;r and independent 

/ m < f - y\ 

- 1 x' J 


i is inversely proportional to 

a h 


- 1 which corresponds to the steady case of Ackeret' s result. In 

Z ° ne 11 (m~+T'= ~"'x' ° =iTT^> varies not with x- or f - t 0 ' 

% - tn 1 

alone, but with the conical coordinate . It is a complicated 


function of M and 


c 0 


Thus the present solution is similar to 


Busemann' s conical flow. So is the behavior of F h = ^ against t ’ - t 


<x h 


0 
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As far 

as 

Oh-™ i« 

a h 

concerned 

in zone 

I. 

In general, 
t* - t 0 ’ 

f h> F h> 

number, 

if 

or 

X' 

f - t 0 ' 


Vert ical Gust 

At time t, if the airfoil begins to encounter a uniform vertical 
gust of velocity g, such -unit-step downwash g at any point on the 
airfoil can also be expressed in terms of Dirichlet' s integral 


v(x,-0,t) = g(t) = g 


sin cn(t - t 0 - g) 


0) 


du> + — 
2 


( 9 ) 


* / ^ \ 

Figure l(b) shows si_J- = a (t) as a function of both t and x. 


U 


Substituting equation (9) into equation (5), the disturbance 
potential 0 g due to vertical' gust can be written as 


0 g (x,-O,t) 


« 2 (m 2 - l) 1 / 2 




r>T, 


•dS 

uq U 0 


2 d r°° sin “(* " t " to " i) 

Jo - T l )( T 2 - T ) 


da) - 


Ox 

uo * u 


PTr 


dT 


T 1 ^( T - T l)( T 2 - t ) 


-g 


jt(m 2 - l) 1 / 2 


dl | ^ sin cult 


> 0 


(t - t Q - | - Mx)j 0 (odX) - 


xg 


2 (m 2 - l) 1 / 2 


(10) 
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or > after integration, 0g has to be also expressed in three zones. 


(m 2 


^[l * -- 1 (5 - >)] 


0 g (x,-O,t) = - t 


(M 2 - l)V 2 


0 


(' 


0 < v < 


M 


(ID 


/ M+l . \ 

\ M = V j 


It is interesting to note that 0h and. 0 g are similar in form. The 
Fourier integral method of Bessel functions used will take care of the 
three zones automatically without considering each zone .individually. 
This is one of the main advantages of the present method. 


Table 3 shows the detailed characteristics of — } -2M. = f } 

a g a g S 

G Lg G Mg 

= — = F g > and - — = G g as functions of time. Figures 5 to 7 show 
a g a g 

fgj F g; and G g as functions of time. The detailed calculations are 
also omitted. 


Pitching 

For convenience , the essential results, of reference 1 for the 
pitching case are given in table 1. Figures 8 to 10 show some charac- 
teristics of this case. 


In the case where xq' = 0 or the airfoil rotates about the leading 

C-nrr 

edge at a constant rate - ag) starting at t' - t 0 ' , f^» = — 

x’ (<*1 " °o) 

. t* - tQ* 

varies linearly with but inversely with M in zone III 


/ t* - tQ 1 • m \ 

\0 = ; = — r ) if the angle of attack is small at any instant. 

\ x 1 M + 1/ 


/ 
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In zone I 


/_m_ s *' - y\ it also 

\M - 1 x' / 


varies linearly with 


t' - t 0 ' 


but 


inversely with - 1. If M > 1.5 f also varies approxi- 

x'(ai - oq ) 


mately linearly with 


t* - t 0 « 


and.' is continuous in value with zones I 


and III although the exact expression as shown in table 1 is quite com- 
plicated. Of course, as the Mach number is near to 1, this linear 
approximation becomes worse. It seems better to use the exact solution 

C lA 

in the table. Similar statements can be madt about Fx = — =r— 

a (*i - &o) 

and G a = -r - x - 

\ a l ~ “o) 

If Xq 1 -f 0 some additional contributions have to be considered as 


fd = 


"pa 


x *(5-i - c^) 


f ‘ 1 - Af • 1 


(12) 


F • = 


-’La, 


= Fa’ - AFa’ 


a ,-r -r \ . a ‘-"■a 

( a l - ao) 


(13) 


G d " 




.a 


(di - do) 


1 = G- 1 - Afi. * 1 

- v 


(lb) 


where 


^d’ 


= f 


AF-* 


h' 


a 


AG- * 
a 


l 0 


F^, and = are given in table 2 and 

V 


figures 2 to L. 
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REACTION ON AN AIRFOIL DUE TO MOTION OF AN ARBITRARY 
TIME -DEPENDENT FUNCTION 


In reference 1 it has been pointed out that, if the solution of a 
linear, ordinary or partial, differential equation with a simple boundary 
condition is known and if the given boundary condition can be obtained 
by superposition of the simple boundary condition, the general solution 
of the same differential equation can then be obtained by superimposing 
the elementary solution corresponding to the simple boundary. In the 
case of the initial -value problem, the statement is also true. With the 
transient problem, the important relation of superposition is the con- 
volution integral. This integral relation is associated with many names 
in modern mathematical physics. In the transient problem of heat transfer 
it is called Duhamel’s integral in honor of the French mathematician ' 

Duhamel (reference 7, pp. to3-lK>4) . In England it is commonly known as • 
Borel s theorem (reference 8, pp. 321-328). In this country it is 
commonly known as Carson's integral, particularly in electrical trans- 
mission. In Germany it is usually called the Faltung theorem (refer- 
ence 9, pp. 159-167), and this was known to Tricomi. The concept of the 
convolution integral plays a very important part in applied mathematics. 

In aerodynamics, Jones applied this integral to the problem of air- 
plane dynamics in incompressible flow (reference 10). The results may 
be considered natural in his problem of ordinary differential equations. 
Garrick (reference 11) applied this integral to find the relation of the 
Wagner Function to the Theodorsen function in the case of two-dimensional 
incompressible nonstationary flow. 


4.4 r ® ference Present author constructed an integral relation 

entirely from the physical concept because in this problem of supersonic 
transient flow the kernel function is discontinuous in the first deriva- 
tive at a certain instant of time. The ordinary concept of the convolu- 
tion integral cannot be applied without detailed examination, although 
the result can be rewritten in quite a similar form to the ordinary 
convolution integral. For example, the pressure coefficient Cp(t',x' ), 

the lift coefficient Cj^t'), and the moment coefficient Cj^(t' ) for 

any arbitrary a(t) = — — — or ~~ can be obtained by means of the 

convolution integral of the kernel function due to unit-step change in 
angle of attack as follows: 


✓ 
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(15) 


where F, G, and f are kernel functions discontinuous in first - 
derivatives as shown in the following table. 


^^Condit ion 
Function^-. 

M St' 
M - 1 

M < t , < M 
M + 1 “ M - 1 

0 % t' < M 

M + 1 

f(f,x’) 

fl(t' ,x' ) 

f 2 (t’ ,x' ) 

f 3 (t’,x' ) 

F(t') 

F]_( t' ) 

F 2 (f) 

F 3 (t' ) 

G(t') 

Gi(t') 

^(t* ) 

G 3 (t' ) 


The subscript i should be equal to 1, 2, or 3 according to the 
above range of time in the term of the left-hand side. Owing to the 
discontinuity of these functions, the integral on the right must be 
broken up into parts according to the range of time for each function. 
For example, in the flapping case, the lift is:. 





(16) 
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(b) For 


M * t , * M 


M + 1 


M - 1 


> ' 


c l/’ t ')ll ~ [o'h^ +0 ^ " %( -oj] F 2 (t' ) + 


U* - 


M 

M+l 


a h ( T ’)F 2 (f - t') dT' + 


n t* 


ft- M 
M+l 


d h (T’)F 3 (t' -• T * ) dr’ 


(17) 


(c) For 0 ^ t* <-i* 

M + 1 ' 


nt' 

C L<V) m = [a^O) - % (-0)]F 3 (f). + a h ( T ’)F 3 (f - t') dr 


( l8 ) 

Tables 2 and 3 give the expressions for f ' s, F’s, and G' s in the 
case of flapping and vertical gusts. In the case' of pitching, the above 
relation is the same as equation ( 15 ), except that ct should be replaced 
by d as given in reference 1 except for C p which is replaced by 

£7 in equation (15). Another complication in this pitching case is the 

location of the axis of rotation Xq’, which is the ratio of the axis 

location from the Reading edge divided by the chord. The contribution 
due to the nonzero value of x 0 * is denoted by Af^, AF^, and AG^ 

and should be added respectively to f & , F^, and G^. As expected, 

Af ^ _ AF^ AGi, 

77T _ f h> ~T = F h> and “ T = Gh as shown in figures 2, 3, and 

0 *0 x 0 

respectively. 

As an application of the above convolution integrals, the cases of 
the harmonically oscillating airfoil starting abruptly from rest at t = 0 
have been investigated as follows: 
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Pitching oscillation about the leading edge (xq’ = 0).- In the case 
of pitching oscillation about the leading edge, ~~~ 


a(t' ) = a for t' < 0 

a(t' ) = ae 1 ^ = ae 1 ^' for t' > 0 


(19) 


where is the frequency of oscillation, p = — , and. t' = — . (See 

u c 

table 4.) Actually, equation (19) is a complex quantity, and the real 
part is 


Real 

a( t ’ ) = a 

for t 1 

< o 

Real 

a(t’ ) = a cos pt' 

for t ! 

> 0 


> ( 20 ) 


The expression a(t' ) = a for t' < 0 is necessary in order that 
d(t') be finite at t' = 0, because only finite d(t' ) is allowed for 
pitching. (For details see reference 1.) The imaginary part is 


Imaginary a(t') = 0 
Imaginary a(t’ ) = d sin pt' 


for t' < 0 
for t' > 0 


( 21 ) 


With the complex a,(t' ), the analysis is more convenient than with 
the real or the imaginary part alone. If the axis of rotation is not 
at the leading i edge Xq’ £ 0, the additional effect can be obtained from 
the flapping case. 

Flapping oscillation .- In the case of flapping oscillation. 


^(t) = 0 


for t < 0 


/ , \ 6 - icut 

a h( t ) = u = a h e 


(22) 


for 


t > 0 
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The solution is given in table 5 for the three time zones. Some of the 
integrals are given in appendix B. 

Harmonically oscillating gust . — In the case of a harmonically 
oscillating gust. 


a g (t) 


a g (t) 


The solution is given in table 6. Some of the integrals are given 
in appendix B. All three cases are essential to the supersonic flutter. 
The present analysis gives the aerodynamic behavior of the airfoil for 
the whole time history, if the harmonic oscillation in flapping and 

pitching start abruptly at t = 0. If t’ > — - — , the present work 

M - 1 

should check with reference 6 exactly. Appendix C shows the comparison. 

To show the present analysis graphically, C L and C M have been 

calculated with M = 1.5 and p = it/2 for the three oscillations. 
Figures 11 and 12 show the and of the pitching oscillation. 

Figures 13 and 14 show them for flapping and figures 15 and 16 show them 
for vertical gust. The corresponding values of a(t’) are also shown 
for comparison of phase shifts. The maximum and are larger 

in the transient beginning than the steady case which is represented by. 

dotted curves. The transient effect dies out completely when t' > M 
• , M — 1 

It is interesting to note that the extremes (maximum or minimum) of 

or Cm always lead the corresponding extremes of a. But in the case 

of flapping oscillation and oscillating gust, the extremes of C L or 

always lag behind the corresponding extremes of a. As the supersonic 
Mach number nears 1, the transient effect becomes more pronounced and 
lasts longer. 

In the time zone III (o S t' ^ ^ , the expressions for C L (t' ) 

and C M (f) are very simple in all three cases of harmonic oscillation 

for pitching, flapping, and vertical gust. (Refer to tables 4, 5, and 6.) 

In the time zone II [ - ** ^ t' ^ — — ), the expressions are rather 

\M + 1 M - 1/ 

complicated and cannot be represented in a closed form. Three new 
functions or integrals have to be defined. Let 


= 0 


t < 0 


(23) 


a g e iaJt , t > 0 
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where m 


M 2 - 1 
M ' 


P > 0, M > 1, 0 = ©i = it , and 6^ 


is usually a 


function of time t' . As designated hy Garrick and Rubinow, J o^ e l) 

may be called the incomplete Bessel function of the zero order. It is 
complex except when, = 0. At 6^^ = 0, equation (2*0 reduces to the 

ordinary Bessel function of the zero order which is real. Simi- 

larly, the incomplete Bessel function of the first order can be defined 
as ' 




cos 6 


exp 


(it cos e ) 


(25) 


which is also complex, but reduces to J-^), the ordinary Bessel 

function of the first order when 6 1 = 0. The new integral is now 
defined as , 


C(p,M; 0i) = 


i T d 0 eJ-^m ) 

n Jq \M - cos 6/ 


(26) 


When 6 ^ = 0^ it reduces to 


C(p 




n * 

f HA £> VDI — 

-ipM A 

J 0 % 

1 - cos 6/ 


(27) 


For convenience, C(p,M) is called the C-function and C(p,M; 0 2 ) 

is called the incomplete C— function. The incomplete C— function occurs 
in time zone II for all three cases. In time zone III, it reduces to 
C(3,M). This function, -to the author's knowledge, has never been 
explored before and is investigated in appendix D. In time zone III, 
Ci^t') and (^(t* ) can be expressed with the known functions 
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except C(f3,M) in all the three cases. It is expected that C(|3,M) • 
is as important to supersonic flutter as the Theodorsen function is to 
flutter in incompressible flow. 


TRANSIENT AERODYNAMIC BEHAVIOR OF THE CONTROL SURFACE 


. Since the principle of superposition holds for the two-dimensional 
linear problem if the deflection angle 5 of the control surface is 
measured from the main wing, the control surface itself may be considered 
as an independent airfoil at the corresponding deflection angle 6 and 
with its chord C s . Under- such a consideration, the lift and moment of 


the control surface itself in nonstationary motion can be obtained, 
directly from the result of the airfoil in the early section, if the 
proper time scale is used. The time required to travel; a length C g 

is t* 


= SSL 


U 


and, if t* is used as the time unit, the nondimens ional 


time t" 

• 1 

time t' 


= — is connected with the true time 

t* c s 

4 = — in the relation 
t C 


and nondimens ional 


(t" - t 0 ")t* = t = t't 


and 



(28) 


where tQ n t* is the difference in starting time of the control movement 
from that of the main wing. Thus for an individual increment in 6 as ' 
the step function at time tj", the lift and moment coefficients of the 
control surface at later times are' 


C L (t") - - V) 

(J = 1,2,3) 

-H") 


(29) 


where Fj and Gj are the same as shown in table 1 except that t' is 
replaced by t". 
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(30) 


Then, the increment of the total wing lift coefficient due to the control 
surface is 
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Similarly, the increment of total moment coefficient due to the 
control surface can he found. 

As an example of the effect of the control surface on the total 
lift, the following calculation is made. 

(1) M = 1.5 a = 1000 ft/sec t 1 = 0.02 sec = 5) 

C = 6 ft C S = 2 ft (k c = 1/3) x Q = 0 x x = k ft 

(2) Operation schedule of the main wing: 



a = a = 0 

f < 0 


a = d 

1 — 1 

P 

V)i 

p 

VII 

0 


d = 0 

V <f 

(3) Operation schedule of the 

control surface: 

6(f) 

= 6(f) = 0 

f < 0 

5(f) 

= 2d 

ti ' 

0 £ f £ -±- 


" 3 

6(f) 

= 0 

V < t , < 2t i’ 
3 =t ■ 3 

5(f) 

= -2d 

— — ^ f S t; • 
3 1 

6(f) 

= 8(f ) = 0 

.f >'t 1 ' 


The distribution of d(f ) and 8(f ) against t' and the geometry 
of the airfoil and control surface are given in figure 17. The contri- 
bution of the control surface to total lift, that is, AC^f ), is shown 

in the lower curve; C L (f ) of the main wing is shown in the dotted 
curve which has been given in reference 1. Since the ratio of' 6(f) 
and a(f ) is given, this curve can be used for any arbitrary &. In 

4df _ 

this curve Ctq = — is used instead of a. 

( M 2 - l) 1 / 2 
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DISCUSSION 


For the airfoil in a flow. of constant supersonic speed, the tran- 
sient effect due to pitching, flapping, or vertical gust will damp out 

Q 

in a time period immediately after the change in angle of attack 

U ~ cl 

ceases. Although the result is obtained from the linearized theory, it 
is expected to be approximately true for the nonlinear theory if the 
angle of attack is reasonably small. After that time period, the lift, 
wave drag, and moment become the same as given by Ackeret in steady two- 
dimensional linearized supersonic flow. The transient effect on C p , 

Cl* and Cj^ becomes more pronounced and lasts longer as the supersonic 

Mach number approaches 1. The same is also true of the transient period 
of the harmonic oscillations in pitching, flapping, or vertical gust 
which start from rest at t = 0. 


In the case of pitching with constant rate. 


"Pa 


■'Let 


X'( a i ■- Oq) ( a X - Oq) 


and 


-Md 


( d i - *0) 


can be approximated satisfactorily with a straight line 


in the time zone II, if M > 1.5. When such an approximation is adopted, 
the case of pitching harmonic oscillation can be evaluated very quickly. 
It is also easy to evaluate any arbitrary motion in pitching. In the 
present analysis, the exact expressions in table 3 were used instead of 
the above approximation because no simple approximation can be obtained 
for either flapping or vertical gust. 


In solving the transient problems with the convolution integral, 
one new function 0(0, M) is discovered. For the time from the abrupt 

start ^t > Cp, and Cjj in the. cases of harmonic oscil- 

lations have to be expressed in terms of C(0,M). It seems of comparable 
importance to supersonic flutter as the Theodorsen function is to flutter 
in incompressible flow. In the transient time zone II, where 

" ~~ ~ t M a new function c(0,M; ©i) , called the incomplete 

C-function, occurs which assumes an importance equal to that of the 
incomplete Bessel functions. More complete calculation of C(0,M) seems 
useful for the analysis of supersonic flutter in order to cover wider 
ranges of Mach number and frequency than those given in reference 6 . 
Appendix D gives some of the properties of C(0,M), C(p,M; © 2 ), and 
the incomplete Bessel functions. . 
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The pressure coefficients 

r» 

•'pg 


C P<i . .... C nh 

~(h - Ao) ln Pitching, J. 


m 


flapping, and 3; — in vertical gust are functions of Mach number and 
. a 8 


a conical parameter 


c 0 


x’ 


only. They are analogous to the behavior 


of Busemann's conical flow. As pointed out in reference 1, these results 

c * n f, PP ^ ied t0 a yaWing in ^ inite win g> if the leading edge is ahead 
of the Mach line. 


The effects of additional degrees of freedom such as the control 
surface or servoflaps can be evaluated with the result of the pitching 
and flapping of the main wing as shown under Transient Aerodynamic- 
Behavior of the Control Surface. With the present basic approach, a 
good aeronautical engineer with ingenuity should be able to solve all 
two-dimensional problems of flutter and any other arbitrary motion. 

Of course in the case of complicated time -dependent functions of angle 
of attack some numerical or graphical integration of the convolution 
integral might be necessary as shown in references 1 and 10. 

For the case of the harmonic oscillation with constant maximum 
amplitude with abrupt start at t = 0 the absolute magnitudes of the 
extremes in C L and C M for pitching and flapping degrees of freedom 

are larger in the transient region — - — - S t S ^ £ than those 

M+1U M-1U 

m the steady region jj < t, particularly as the supersonic Mach 

number nears 1. In addition, the maximum and minimum of the oscillating 
angle of attack are not in phase with the maximum and min imum of 

or C M . In the case of’ pitching, the load leads the angle of attack; 
in the case of flapping, the angle of attack leads' the load. Also, the 
angle of attack leads the load in the case of vertical gust. 

As an interesting example of the harmonic oscillation building up 
to flutter or damping out, the case with complex should be inves- 
tigated. With the convolution integral, it seems within the reach of 
the present analysis. Of course, the imaginary part of cu must be 
determined from the interaction of aerodynamic forces and the elastic 
behavior of the wing. 


The Johns Hopkins University 

Baltimore, Md. , November 4, 1949 
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APPENDIX A 
SYMBOLS 


a 


C 


C LO 


-M 


C S 


velocity of sound 
wing chord 

transient lift coefficient 

■4 

two-dimensional lift coefficient in steady case 
transient moment coefficient, about axis of rotation 
pressure coefficient ^2(p - 
chord of control surface 



C(P,M) = 




d0 exp 


-ip 


1 - 


cos 6 
M < 


= i f* de exp (— -m ) 

* Jo V M " cos Q ' 



F( ) 
f( ) 
G( ) 

g 

g 


kernel function concerning lift coefficient 
kernel function concerning pressure coefficient 
kernel function concerning moment coefficient 
velocity of vertical gust 
•maximum velocity of uniform vertical gust 



22 


NACA TN 2333 


fi descending velocity 

fi maximum velocity of vertical flapping wing 

J 0 ( ) Bessel function of zero order 

J^( ) Bessel function of first order 



M Mach number 


m 


parameter of Mach number 



t 

t = c/u 


t* 


t* = Cg/u 

t" 

U 


v 

X 

X’ 

*0 

y 


a 

a 

v 5/u 

a h = fi/u 


time 


nondimens ional time 



nondimens ional time (control surface) (t/t*) 

free-stream velocity 

y-component of velocity 

axis along chord direction 

percent of chord (x/c) 

axis of rotation of entire airfoil 

vertical axis 

angle of attack 

constant angle of attack 


/ 
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3 frequency parameter for oscillating airfoils 

5 deflection angle of control surface 

t] source location along y-axis 

6 running variable 
lower limit of 6 

X = x/a(M 2 - l) 


v - x/u(t - t 0 )= x*/(t* - t Q ’) 

£ source location along x-axis 


a = J3M 

T 


0 


CD 

Subscripts: 


time interval 

nondimens ional time interval (Vt) 
velocity potential 

frequency parameter (angular velocity per second) 


g 

h 

d 


vertical gust 
flapping 

pitching (changing angle of attack) 


5 


control surface 
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APPENDIX B 
SUMMARY OF INTEGRALS 

Important Integrals in the Case of Steady Harmonic Oscillation 

The following integrals are used in evaluating C L and C M for the 
time zone t* > - M (as shown in tables k, 5 , and 6); 
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PIS 
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Important Integrals in the Case of Transient 
Harmonic Oscillation 

• I 

The following integrals are used in calculating C L and C M for 

the time zone — M ^ t’ ^ — M (as shown in tables 4, 5, and 6). 

M + ± M - I 
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M+1 e^ T, (t' - t 1 ) 3 cos" 1 


t ,:5 , M( t 1 - 1) 

cos x 

IP t 1 


(M 2 - 1 ) 


M 3 e P ' 

•n 5/2 « 


!*■-".) r 


os "* 1 ( mt’ -M) 


^M 2 + 1 + 2M cos 0 + ^ cos 0^ 


£OS 0 

e“ iP T~ d0 + 


— 7- \|l - (mt' - M) 2 + — I 7 ' 

W („2 . i)3/2 V iP 7 


*•- » 


M * 1 eiP T '(t' - t') 2 \Z(t' - 


- T ' -1)2 dr' - *£ jl - (mV - M) 2 f 

(M 2 - 1 ) 5 /? ^ L 


— (M + 3n>t' ) + 6(±fl ^ V ~_^ I 

1P (M2 -1)5/2 l P J cos -l (mtt . M) 


r 

• 2M + ^ + cos 0 [m 2 + 1 + — 

-ifmt.-M! 13 L 13 


2 l -13 


..cos^ 
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APPENDIX C 

CHECK WITH NACA TN 1158 ' 


To check the present results through convolution integrals with 
Garrick and Rubinow 1 s work, a case of the steady harmonic oscillation 
of the angle of attack about the leading edge is investigated. With 
their notation, C^ due to the harmonic oscillation of the angle of 

attack alone /a = a^e may be verified from their equation ( 26 ) as 


C KX ' + iL h) < C1 > 

At x Q = 0, 

= l 3 = ^3* “ 2 xqLi 

Lj| = Lj^' Lj^ = Ljj_' - 2 xqL2 


At M = 2, for 


— = 1.661, from their table II, 
k ’ 


l 3 ’ = 1.50219 L 4 ' = 0.69968 


Hence, 


icat 


c La = .ao(2.l6315 + 1.0075^i)e 
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NOTATION 


NACA TN 1158 

Present paper 

t 

t’=™ 


C 

0) 

— 

ZD 

II 


U 

b 

C = 2b 

V 

> 

II 

II 

!z> 

k = ^ 

p = 2k 

V 


icut 

e 

exp (ipt') = exp (iojt) 


In the present case a = a exp (ipt' ). From table 4 for t' ^ M 

M - l' 


r 


— — — = ip exp 

Llv . * Y 


jM 

' -1 
r— 1 

i 

OJ 

■ -♦ 
OJ 

1 

a 

L a « 2 - J 


0\mJ 

Jm 2 -1 M^M 2 - 1 iP J 



M 3 ( ip) 2 


[l - C( P,M)j 


whence/ at M = 2 and p = 1.2, C L = a exp ( ipt' >(2.16315 + 1.00753i). 
The expressions are identical, within computational limits. 
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APPENDIX D 

INVESTIGATION OF THE FEW IMPORTANT INTEGRALS RELATED TO 
SUPERSONIC FLUTTER AND TRANSIENT PROBLEMS 
C -Function 


The problem is investigation of the C-function or the integral 


C(P,M) = 


if™ / -iPM \ 

= - ae exp 

« \m - cos e 

uo ' ' 


(P > o) 

(M > l) 


If the above integral is differentiated with respect to 3 and M, 


dC 

dM 




cos 6 


(M - cos e)‘ 


(Dl) 


= — I ae [exp 


dc _ -i 

dp * 


L 


/ -1PM \ 
\M - cos ej 


M 


M - cos e 


(D2) 


d 2 C 

dp 2 


T r *[~(r^ 


w 1 


(M - cos e)‘ 


(D3) 


With the above three relations, it is found that C(p,M) 
satisfies the differential equation 


- Co + iCn 


dC = i cTC + M dC 

an " p 3 m 


(dm 



32 


MCA TN 2333 


If the real and imaginary parts of C(3,M) are taken separately, 

d g R _ _ + M SCr 

Sp " ^2 3 bu 


^1 _ , M dC T 

dp dp 2 3 & 


(95) 


which are two simultaneous differential equations. The required six 
boundary conditions can be obtained from C(p,M) as follows. 

For 3=0, 


For M — > oo. 


c r (o,m) = 1 

Cp( o,m) 

bc R _ 

1 

H 

O 

ro 

dp 

as " f 

= cos P 



-M 


(D6) 


(D7) 


The integral ± J d0 exp 


/ -1PM \ 
\M - cos 3/ 


can be written as 


i f” de exp l— l ° ) 

n Jo \ M - cos e ) 


C(a,M) = ~ I <10 exp 

Jo 

where a = PM. Under the new definition, the function C!(a,M) can be 
shown to satisfy the differential equation 


ia 


d 2 C _ be 

bo 2 


(D8) 
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or . . 

o = 

do 2 _8M 

> 

h 2 C r dC R 

a i = - — — 

So 2 

It is easy to shov that the real and imaginary parts of 
must satisfy the equations 


and 



The corresponding boundary conditions are: 
For a = 0, M > 1, 


Cr = 1 

dCt 


dcr 


= 0 


d 2 Cj 

ha 2 


C T = 0 


he 

ha 


I _ -1 


- 1 


M 


(m 2 - 1) 3 / 2 


S 2 C] 


= 0 


For M — > co, 


txi 

il 

H 

Ci = 

dc R 

^ = 0 
dM 

SCj 

dM~ 

dC R 

h^~ = ° 

dCj 



> 


(D9) 


C(a,M) 


(DIO) 


(Dll) 


(D12) 


0 



The preceding function is so important to the supersonic flutter 
problem that thorough investigation of the function and the other 
solutions of the corresponding differential equation seems advisable. 

To show the general behavior of this new function, a case with M = 1. 3 
and 0^3 = 3 has been calculated numerically. See table 7 and 
figure 18 . 


Incomplete C-f unction 
The incomplete C-function is 


C(P,M; 


0 l) 



exp 


\M - cos 6 ) 


de 


(D13) 


It can be shown that this incomplete C-function also satisfies the 
differential equations (D5) except the boundary conditions are obtained 
from equation (D13). 

For 3 = 0, 




For M — > oo. 



> (D14) 



J 





COS 3 



(D150 


To show the nature of the incomplete C-function, a case with M = 1. 
3 = ana 0 ^ ^ « has been calculated by means -of numerical inte- 

gration. See table 8 and figure 19. The above data are used in the tran 
sient behavior as shown in figures 11 to 16 . 
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Incomplete Bessel Functions 
The incomplete Bessel functions are 


On 

= — exp (+iz cos 6) d 6 (Dl6) 

J e l 


cos 0 exp (+iz cos 6) d G (D17) 

These two functions become Jq(z) and J-^(z), respectively, when 
= 0. They have been called the incomplete Bessel functions by 
Garrick and Rubinow. The case M = 1.5* P = and 0 ^ ^ jt has 

been calculated by numerical methods. It is interesting that each has 
a real and an imaginary part. At 0^ = 0, both become real and ordinary 
Bessel functions. Both play as important roles as the incomplete 
C-function in harmonic oscillations at the transient time zone II. See 
table 8 and figures 20 and' 21. 
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CHARACTERISTICS OF A FLAPPING AIRFOIL WITH li (t > 0) OR ADDITIONAL 
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CHARACTERISTICS OF AH AIRFOIL ENCOUNTERING A SHARP VERTICAL 
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TABLE 




CHANGING ANGLE OF ATTACK 


E l - ae 1 ^’ t’ > 0 

in the caBe Xq' =* 0. See flge. 11 and 12. 

l - a ' 0 > t* 



























VERTICAL GUST 
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C-function. 

Tncomplate C-function C(p,M; 0^; see appendix D. 
Incomplete Bessel functions; see appendix D. 
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.3 


3 

Heal 

Imaginary 

0 

1.0 

0 

.3 

• .84153 

-.4l796i 

.5 

. 62101 

-.56797i 

.8 ' 

. 32000 

-.569621 

1.0 

.21844 

_.49844i 

1.2 

.19147 

-.44854i 

1.5 

.17883 

-.477121 

1.7 

.12058 

-.534231 

,2.0 

-.48081 

-.562791 

2.2 

-.16005 

,-.511541 

2.4 

-.22697 

-.4249li 

2.5 

-.23982 

-.381231 

2.93 

-.22265 

-.293551 
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TABLE 8 

VALUES OF INCOMPLETE C -FUNCTION AND INCOMPLETE BESSEL FUNCTIONS 


. [p = '«/2, ■ M = 1.5] 


«1 

c(p,M 

; e i) 

J o(i 

r> e i) 


0 

-0.05339 

,-0.47l48i 

0.29055 

0 

0.58148 

0 

.2 

-.04952 

• -.534931 

. 30946 

.060781 

.52110 

. 018801 

.4 

-.02526 

-.593171 

.32378 

.122781 

. 46198 

.032491 

.6 

. 02500 

-.630791 

. 32891 

. 1861 4i 

.40647 

.037071 

.8 

.08770 

-.634931 

. 32058 

.249101 

. 35837 

.030881 

1.0 

.14434 

-.60720i 

.29582 

.307501 

. 32204 

.015731 

1.2 • 

.18435 

-.558191 

.25417 

.355361 

. 30020 

-.002891 

1.4 

.20575 

498 *18 i J 

.19871 

. 385801 

.29172 

-.015721 

1.6 

.21084 

-.435171 

.13598 

. 394201 

.29074 

•r. 021901 

1.8 

.20310 

-.372061 

. 07454 

.379001 

.28841 

-.014091 

2.0 

. 18569 

-.310871 

.02257 

. 3428li 

.27643 

' .002481 

2.2 

.16124 

-.252121 

-.01429 

.291251 

.25027 

.020821 

2.4 

.13173 

-.195731 

-.03407 

.230961 

.21010 

.033771 

2.6 

.09867 

-.141331 

-.03808 

.167561 . 

.15938 

.036831 

2.8 

.06326 

-,08844i 

-.02980 

. 104501 

.10248. 

.029291 

3.0 

.02644 

-.036501 

-.01367 

.042951 

.04279 

.013591 

3.1 

.00778 

-.0107H 

-.oo4o8 

.012591 

. 01259 

.004111 

1 

Tt i 

0 

0 

0 

0 

0 

1 

0 

1 
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Figure 8.“ Plot of f F for xF - 0 as a function of - — at various Mach numbers. 








Figure 9.- Plot of Fi' for Xp,' = 0 as a function of t' - t n ' at various Mach numbers. 
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t' 


Figure 17.- Distribution of d (t 1 ) and 6(t') against t' and geometry of 
airfoil and control surface. Contribution of control surface to total lift 
shown by curve at bottom of figure. CjJt') shown by dotted curve 
(from reference 1). tj_' = 5. 
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Figure 19.- Incomplete C -function C(p,M; 0 ^) for 
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Figure 21.- Incomplete Bessel function of first order for 

M = 1.5 and p = n/2. 
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